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We discuss the problem of universe acceleration driven by global rotation. The redshift-magnitude 
relation is calculated and discussed in the context of SN la observation data. It is shown that 
the dynamics of considered problem is equivalent to the Friedmann model with additional non- 
interacting fluid with negative pressure. We demonstrate that the universe acceleration increase is 
due to the presence of global rotation effects, although the cosmological constant is still required 
to explain the SN la data. We discuss some observational constraints coming from SN la imposed 
on the behaviour of the homogeneous Newtonian universe in which matter rotates relative local 
gyroscopes. In the Newtonian theory fi r ,o can be identified with W] o (only dust fluid is admissible) 
and rotation can exist with f2 r ,o = fi u ,o < 0. However, the best-fit flat model is the model without 
rotation, i.e., fi Wj o = 0. In the considered case we obtain the limit for fi^o > —0.033 at la level. 
We are also beyond the model and postulate the existence of additional matter which scales like 
radiation matter and then analyse how that model fits the SN la data. In this case the limits 
on rotation coming from BBN and CMB anisotropies are also obtained. If we assume that the 
current estimates are f2 m ,o ~ 0.3, f2 r ,o ~ 10 -4 , then the SN la data show that fi w ,o > —0.01 (or 
ujq < 2.6 ■ 10~~ 19 rad/s). The statistical analysis gives us that the interval for any matter scaling like 
radiation is n r , G (-0.01,0.04). 



I. INTRODUCTION 



We consider a homogeneous universe which is more general than the FRW model, in which matter additionally 
rotates relative to localgyroscopes pj . The motion of the fluid in such a universe is described by the scalar expansion 
9, the rotation tensor u a b, and the shear tensor a a b- The homogeneous rotation of fluid as a whole is usually called 
the global rotation of the universe 0] . 

Applying these concepts we must remember that CMB strongly restricts (indirectly from observations) the value 
of angular velocity 0, On the other hand one can estimate the present value of 9$, ^0: and ctq also directly from 
observations of galaxies 5j. These observations show that 6* = 3H , ui n < %/%_ 0o ~ #o/4- 

The propagation equation for 9, known as the Raychaudhuri equation [a, 0], for the perfect fluid with energy- 
momentum tensor T a b = (p + p)u a Ub + pg a b (where p and p are the energy density and pressure, respectively), is 



9 - «« + Je 2 + 2(a 2 - u 2 ) + ~(p + 3p) - A = (1) 

where ii a = u a -,bU b is the acceleration vector; we shall use a dot to denote the rate of change of any quantity as 
measured by an observer moving with 4-velocity u a ; and co 2 = aj a t,uj ab /2, a 2 = <J a b<J ah /2 are the scalars of rotation 
and shear, respectively; and A is the cosmological constant. 

If we define a representative length I along a particle world line by 0, Q 

then / represents the volume behaviour of the fluid completely. For example from I one can define the Hubble function 
H and the deceleration parameter q by 
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Using definition (J2J) and @, equation can be rewritten in the form 

3j = 2(u 2 -o a )+u a . a -±(p + 3p) + A. (4) 

This shows how the acceleration of the universe (the curvature of curve l(t)) is directly determined at each point of 
spacetime. Let us note that A acts as a constant repulsive force whereas rotation as a variable repulsive force. 

When oj 2 , a 2 , and u a . a are given as a function of I we can integrate equation (@J. To simplify matter we take u a — 
(because the acceleration vector represents the effects of non-gravitational forces it vanishes when a particle moves 
along a geodesic, which would necessarily follow in the case of dust) . 

It has been shown that spatially homogeneous, rotating, and expanding universes with the perfect fluid have the 
non- vanishing shear This is quite contrary to the case of the homogeneous Newtonian cosmology where many 

such solutions are known. These homogeneous shear-free solutions are independent of the pressure which may be 
set equal to zero or a constant. This difference in the two theories seems to be both surprising and interesting since 
Ellis' theorem has a purely local character, and it is completely independent of the strength of the gravitational field 



Let us consider solutions with a = = ii, ujQ = 0. In this case uj 2 — uj 2 /1 4 where uj — 0. Then we can integrate 
the Raychaudhuri equation using the conservation equation 

P + Q(p + P ) = 0. (5) 



The occurrence of term p in the factor {p + p) is a special relativistic effect 

In the considered case of wO ^ we obtain the generalised Friedmann equation 

3/2 _|_ _ t^. _ \p = _ 3fe (6a) 
l z I 

k = (6b) 

where p = const, p = and l(t) = a(t) is the scale factor. From the mathematical point of view equation © is a first 
integral of system |@} ■ 

Equation © can be treated as basic equations in a Newtonian homogeneous cosmology. Solutions of this equation 
represent shear-free Newtonian cosmologies which are in general both expanding (i ^= 0) and rotating (u)q ^ 0). 
Equation |J?J| is called the Heckmann-Schiiking equation 0. 

If we consider models with w = = u, cr9 ^ 0, and the Ricci tensor 3 R a b is isotropic then we obtain a 2 = S 2 // 6 
where S = 0. We can then integrate the Raychaudhuri equation to obtain the generalised Friedmann equation 

Y 2 ill 3 

3i2 ~l±~l~ Af = " 3fc (7a) 
k = (7b) 

where I = (aia2as) 1 / 3 is an average scale factor. 

Therefore, it seems reasonable to assume that a is sufficiently small compared with u> since the shear falls off more 

rapidly than the rotation |2j, |6J, [jj . 

In the case of dust a 2 oc a~ 6 whereas oj 2 oc a~ 4 . The conservation of angular momentum gives u>pa 5 = const 
From equation (J7J we see that the effect of anisotropy is like in the FRW model with stiff matter. In our further 

analysis of observational effects we consider equation © as a simplest model in which the effect of global rotation 

can be investigated. However, we also consider the presence of additional non-interacting radiation matter which can 

be treated as a simple extension beyond the Newtonian model. 



II. EFFECT OF GLOBAL ROTATION ON ACCELERATION OF THE UNIVERSE 

The supernovae observations indicate that the Universe's expansion has started to accelerate during recent cos- 
mological times, and CMB observations su gge st that the Universe is dominated by a dark energy component, with 
negative pressure, driving the acceleration |lCt ITU . While the most obvious candidate for such a component is the 
vacuum energy a plausible alternative is the dynamical vacuum energy or quintessence. However, these models usually 
face fine-tuning problems, because there is a question of explaining why the vacuum energy dominates the Universe 
only recently |12|. 
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To study the effect of the global rotation on the acceleration of the universe we formally introduce rotation to the 
model by definition of 

Pcfr = Pm + Puj = Pmoa~ 3 + Pmo-~ a + A (8a) 
Pcff = ^Pu - A (8b) 

where p w0 — — 2wq < and p w = (like for radiation matter). 

Therefore, in the case of dust filled universe, the dynamical effect of global rotation is equivalent to an additional 
non-interacting fluid with negative pressure. 

In order to take into account the effects of rotation we introduce 

n w= ^_ = _^|f^y 4 (9a) 

m 2 m 2 \a ) y ' 

n Pm POm / a \ 

n - = 3^ = 3HfUJ • (9b) 
For our purpose it is also useful to rewrite the dynamical equations to a new form using dimensionless quantities 

x=—, T=\H \t 
a 

with H = a/a, p cr ,o = 3i?Q and the subscript means that a quantity with this subscript is evaluated today (at time 
to). Additionally we define Qj^o = —ik/QH^ and fl\ t o = A/3-ffp. 
The basic dynamical equations are then rewritten as 

i 

x = -^Y, n ^(l + 3 1 i)x- 2 ~ 3 ~<> (10b) 

i 

where i = (m,w, A). The above equations can be represented as the two-dimensional dynamical system 

x = y (11a) 

y = -lj2 n ^ 1 + 3 ^ x ~ 2 ~ 3Wz ( llb ) 

i 

or by the Hamiltonian dynamical system with the Hamiltonian given in the form 

H = ^x 2 + V(x) = (12) 

and with the potential 

v(x) = -in fc , - tkox- 1 -*". 

i 

The system should be considered on the zero energy level. 

The form of l|12|l can be useful in particle-like description for the simplest model with global rotation, whereas form 
(|llfl is helpful in the analysis of dynamics on a phase plane (x,y). 

The system under consideration can be identified after taking 

Wi = 1/3 (effect of rotation or radiation) 

W2 = (effect of dust matter) 

W3 = — 1 (effect of A, p = A) 

As an example of application of these equations consider the case of O mi o, ^ w ,0; ^A,o 7^ 0. Then our Universe 
accelerates provided that the potential V is a decreasing function of its argument 

= -flu.ox^ 3 - ^O m oar 2 + Qa,qx > 0, (13) 
ax 2 
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i.e., if il m .o = 0, the universe always accelerates for every x. 
For Oa o = the Universe accelerates provided that 

x < -- — 

"m,0 

where £1^,0 < and J2i ^i,o + ^fc,o = 1 (i = m,w, A). 
Our Universe accelerates at present provided that 

-2 - il m .a + 4S1 A>0 + 20 fe ,o > 

or £l\,o > 0.425 where we assume £lk,o — 0, £l m ,o = 0.3. We can see that rotation lowers the value of cosmological 
constant needed to explain the SN la data. 



III. MAGNITUDE- REDSHIFT RELATION IN THE MODEL 



The important test to verify whether rotation may represent "dark energy" (which can be called true dark radiation 
because causes the acceleration of the Universe) is to compare rotation effects with the supernovae type la data. The 
answer is that global rotation may be seriously taken as a candidate to describe only part of dark energy and the 
cosmological constant is still required. 

It is well known that cosmic distance measures, like the luminosity distance, depend sensitively on the spatial 
geometry (curvature) and dynamics. Therefore, luminosity depends on the present densities of the different com- 
ponents of matter content and their equations of state. For this reason, the magnitude-redshift relation for distant 
objects is proposed as a potential test for cosmological models and play an important role in determining cosmological 
parameters. 

Let us consider an observer located at r = at the moment t = t Q who receives light emitted at t = t\ from a 
source of absolute luminosity L located at the radial distance r\ . Of course the cosmological redshift z of the source 
is related with t\ and to by the relation 1 + z = a(to) j a{t\) ■ If the apparent luminosity of the source measured by the 
observer is I, the luminosity distance of the source, defined by 

1 = A (m) 

is 

d L = (1 + z)a n. (15) 

For historical reasons, the observed and absolute luminosities are defined in terms of K-corrected observed and 
absolute magnitudes m and M, respectively (I = 10~ 2m / 5 ■ 2.52 • 10~ 5 ergem" 2 s~ 2 , L = 1(T 2M / 5 • 3.02 • 10 35 ergs" 2 ) 
[l3|. When written in terms of m and M, equation (|14H yields 

m(z,M, n m . ,nA,o) =M + 5log 10 [V L (z,n m! o,n A>0 )] (16) 

where 

M = M -5 log 10 H + 25 (17) 

and 

T> L (z,Q mfi ,Q A:0 ) = H Q d L (z,Q mfi ,Q Afi , H a ) 

is the dimensionless luminosity distance in Mpc. 

By using expression for the FRW spacetime metric we obtain coordinate distance r\ , appearing in (|15H 

Hn)= J d 4 = -j^=, (18) 

J aa J VI - kr 2 

a /(l+z) r 
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with 

tp(ri) = sin -1 r\ for k = +1 

V»(ri) = ri for fe = (19) 

V'( r i) = sinh -1 7* i for fe = — 1. 

By using Hamiltonian constraint l |12(l for the model with dust matter, cosmological constant, curvature, and global 
rotation we obtain 

z 

ip{rx) = I [n fc , (l + ^') 2 + On,o(l + -z') 3 + ^.o(l + z'f + riA.o}- 1 / 2 ^'. (20) 



We obtain finally 

v L ({z, n m , , n A , , n w>0 ) = I V£ / [(i - ^m,o - a,,o - n A ,o)(i + *'Y 



/c 



where 



+ n m>0 (i + z'f + n Ui0 (i + z'f + ft A , r 1/2 ^' , (2i) 



£(x) — sin a; with /C = — f2fe,o when fi^o < 

£(x) = x with JC = 1 when Sl^.o = (22) 

£(x) = sinhx with JC = Qk,o when Qk.o > 



and 



k 

Thus for given A4, fi mi o, ^a,0; S^fc.o, ^,0i equations (|16[1 and (|21|) give the predicted value of m(z) at a given 2. 



IV. MAGNITUDE- REDSHIFT RELATION IN THE MODEL — RESULTS 



We decided to test our model using the Perlmutter sample [10| . To avoid any possible selection effect we choose the 
full sample without excluding any supernova from that sample. It means that our basic sample is Perlmutter sample 
A. We test our model using the likelihood method [rl|. 

Firstly, we should estimate the value of M. (equation Q17))) from the full sample of 60 supernovae taking f^.o = 
(the pure Perlmutter & Riess model). We obtain value of M, = —3.39 (we also assume that the present value of the 
Hubble constant is Hq — 65kms~ 1 Mpc -1 ) what is in very good agreement with the result of Efstathiou et al. 
and Vishwakarma 15) (Vishwakarma obtains A4 C = 24.03 for c = 1, i.e., M. = —3.365). Also the value of \ 2 obtain 
for Perlmutter's flat model is 96.5 what is in very good agreement with Perlmutter's results (see Table 3 in fiojp. 
(Some marginal differences are probably because in our analysis we include both errors in measurements of magnitude 
and radial distances). 

We consider the pure Newtonian model with J\,.o < and assume that fl m ,o ~ 0.3 [l6l . ll7fl . Using the minimalization 
procedure, described below, with aforementioned assumptions we obtain the density distribution for f2 Wj o- The results 
are presented on Fig. ^ Here we find that the limit for fl^.o > —0.033 on the la level, while fl u ,o > —0.065 on the 
2a level. 

The analysis of the pure Newtonian model is presented on Fig. [3 with the magnitude-redshift relation for real data 
(marked with asterisks) and for predicted values by models. The top line is the pure Perlmutter flat model with 
^m,o = 0.28, f2 Aj o = 0.72. The bottom line is the pure flat model with the cosmological constant fl\,o = 0. Between 
these models there are located our models with Q u .a = —0.01 best-fitted model (lower curve) and best-fitted flat model 
(upper curve). The latter model curve overlaps the Perlmutter model curve. One could observe that the difference 
between our lower best-fitted model and the Einstein-de Sitter model with f^A.o — is the largest for z between 
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0. 03 




FIG. 1: The density distribution for Q w ,o in the Newtonian model. We obtain the limit fi W] o > —0.03 on the la level, while 
fiw.o > —0.06 on the la level. 

0.6 and 0.7 and significantly decreases for higher redshifts. There are significant differences between predictions of 
these models and Perlmutter's one where differences to the pure flat model increase for higher redshifts. It gives us 
possibility to discriminate between the Perlmutter model and our model when data from supernovae more distant 
than z ~ 1 could be available. It is very important because for present data our model is only marginally better than 
the Perlmutter model. 

We can also admit that the total matter content scales like radiation. It means that the contribution coming from 
fi^o is included in O r; o- Therefore, in the more detailed analysis we assumed that Qk.o G [ — 1, 1], f2 m ,o G [0, 1] . 
From the formal point of view then we obtain the best fit (x 2 = 94.7) for f^o = —1-0, £l m ,o = 0.54, fl^.o = 0.15, 
^A,o = 1-31, which is completely unrealistic. However, we should note that we obtain, in fact, a three-dimensional 
ellipsoid of possible models depending on O mi o, fi w ,0j ^A,o- It is more complicated than in the case of Perlmutter's 
analysis when he obtains only two-dimensional ellipsoid (depends only on f2 m .o and ^a,o)- But, knowing the best-fit 
values has no enough scientific relevance, if not also confidence levels for parameter intervals are presented. On the 
Fig.[3]we show the levels of constant \ 2 on the plane (^a,0j ^m,o) minimalized over the rest of parameters. The figure 
shows the preferred value of f2A,Oj^m,o- The minimalization procedure confirms the chosen value of M. — —3.39, 
because it is a best-fitted value for the flat models. Now, we would like to obtain confidence contours in the Oa,o> ^m,o 
plane. 

Since from the formal point of view we have no a priori constraints on cosmological parameters we assume here 
that fifc o an d ^m,o are of any value. The result of our analysis are presented on the Fig. This figure shows the 
confidence levels of 2-dimensional distribution of (f2 m ,o, ^A,o)- It is analogous to the confidence level figure obtained 
by Perlmutter. 

Another considered case is the flat model (O^o = 0) where we obtain "corridors" of possible models (we presented 
confidence contours in plane Fig.0. The formal best-fitted flat model is O m ,o = 0.12, £1^ = 0.12, flA,o = 0.76, 
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FIG. 2: Residuals between the Enstein-de Sitter model and three cases: the Einstein-de Sitter itself (zero line) the best-fitted 
model with f2 w> o = — 0.01 (middle curve), and the best-fitted flat model with Q w> o = —0.01 and the Perlmutter model (these 
two curves overlap) (highest curve). Therefore in the case of dust matter, the difference between the Perlmutter model and 
best fit model with rotation becomes detectable for redshifts z > 1.2. 



X 2 = 95.7. In probably a more realistic case we obtain for flat model fi m) o = 0.28, flu>,o = 0.02, i.e., Qa.o = 0.73. For 
that model x 2 = 95.9. For the flat model with low rotation ft mj0 = 0.33,' fi w>0 = -0.01, i.e., tt A . = 0.68, x 2 = 96.0. 
The value of x 2 is practically the same in all three cases. It clearly shows that statistical analysis is not sufficient 
for discrimination between statistically available models. To choose the physically plausible model we need aditional 
information which can be obtained for example from extragalactic astronomy investigations (especially estimations 
for O mi o and Slk,o are useful). 

It is interesting to observe how presence of non-zero (but rather small) Qkfi with realistic rotation changes our 
results. For example for fife o = 0.1, fi m o = 0.28, H u q = —0.01, i.e. CIa,o = 0.63 x 2 = 96.2, while for Sl^ o = —0.1, 
fim.o = 0.38, fi w o — —0.01, i.e. f^A.o = 0.63, x 2 = 95.8. It shows interesting possibility if separately we could find 
value for rotation fi^.o and matter n mj o, than we could test the value of fi* o more precisely than with the models 
without rotation. 

From Ref. [PoIITtT] we obtain that value fl m .o should be not far from 0.3. With this assumption we could find from 
Fig. that fi^o should satisfy fi W) o > —0.01 which gives critical angular velocity ujq = 2.6 • 10 _19 rad/s, is in a good 
agreement with other limits, however it should be pointed that our limit is weaker, Li suggested loq = 6- 10 -21 rad/s. 
In terms of density parameter that limit requires $7^,0 > —5.3 • 10 -6 whereas to obtain Ciufolini and Wheeler's limit 
is required to U)0 > -1-4 ■ 10~ 4 . 

One should note that we give our analysis without excluding any supernovae from Pcrlmutter's data. However, 
from formal point of view, when we analyse full Pcrlmutter's sample A, all analysed models should be rejected even 
on the confidence level 0.99. One of the reasons could be the fact that assumed errors of measurements are too low. 
Nevertheless, another solution is usually suggested. We can exclude 2 supernovae as outliers and 2 as likely reddened 
ones from the sample of 42 high redshift supernovae and eventually 2 outliers from the sample of 18 low redshift 
supernovae (Perlmutters's sample B and C, respectively). We decided to use full Perlmuttler's sample A as our basic 
sample because rejecting any supernovae from the sample could be the source of not fully controlled selection effect. 
On the other side such procedure also could be useful. It is the reason that we decided to check our analysis using 
Perlmutter's samples B and C. It does not significantly change our result, but increases quality of the fit. The formal 
best-fit for sample B (56 supernovae) is (x 2 = 57.5) what gives Qkfi = —0.3, fi m ,o = 0.2, f7 W:0 = 0.17, i.e., Qa,o = 0.93. 
For the flat model we obtain (x 2 = 57.6) £l m .o = 0.03, Q W) o = 0.19, i.e., 17a, o = 0.78, while for "realistic" model 
(Q mfi = 0.28, Q u<0 = 0.03) A ,o = 0.69 x 2 = 57.7. For the flat model with small rotation Q Ut0 = -0.01, ft m . = 0.34, 
i.e., fi A ,o = 0.67, x 2 = 57.8. 
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FIG. 3: Levels of constant \ 2 on the plane (fi m ,o, £)a,o) minimalized over the rest of parameters. The figure shows the preferred 
value of fiA,o, £2 m ,o- 



The formal best-fit for sample C (54 supernovae) (v 2 = 53.6) gives Slk,o = —0.1, Sl m ,o = 0.11, fi W) o = 0.18, i.e., 
Oao = 0.81, while for flat model f2 m ,o = 0.05, O w o = 0.19, i.e., fl\o = 0.86, \ 2 = 53.6, while for "realistic" model 

2 



Q m , = 0.24, fi^.o = 0.07, i.e., fi A ,o = 0.84, x = 53.6. For the flat model with small rotation Q^o = -0.01, 
O m ^ = 0.36, i.e., h A ,o = 0.65, x 2 = 53.7. 

It again confirms our conclusion that on the base of pure statistical analysis we could only select "corridor" of 
possible models. However, if we assume that the Universe is flat Qk.o = 0, we obtain estimations for f2 m .o, S\,,o what 
seems to be realistic. 

One should note that we also could separately estimate the value of M. for sample B and C. We obtain M. = —3.42 
what is again in very good agreement with result of Efstathiou et al. |14| (what for the "combined" sample obtain 
the value of M — —3.45). However, if we use that value in our analysis it does not change significantly our results 
(value of x 2 does not change more then 1 what is marginal effect for x 2 distribution for 53 or 55 degrees of freedom. 

We also analyse the influence of rotation for the age of the Universe. The results are presented on Fig. HJ1 If we 
assumed that f2 m ,o = 0.3 and Hq = 65km/sMpc then small rotation fi^o = —0.01 increases the age of the Universe 
from 14.57 • 10 10 yr to 15.17 • 10 10 yr. 
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FIG. 4: Confidence levels on the plane (f2 m ,o, Ha.o) minimalized over the rest of parameters for the flat model. The figure 
shows the ellipsoid of the preferred value of f2 m ,o, SX\,o- 



Finally, let us study the angular diameter test for our universe. The angular diameter of a galaxy is defined by 



(23) 



where d is a linear size of the galaxy. In a pure flat dust model universe has the minimum value for z m i n = 5/4. It 
is particularly interesting to notice that for flat models with Ha,o 7^ the dark radiation can increase the minimum 
value of 9 toward the largest z m j n and smaller m i n , while the ordinary radiation lowers this value. 

We presents influence of rotation for the angular diameter 0(z) as a function of redshift z. For the flat model with 
^m.o — 0.3 as shown in Fig. [7| the rotation causes the minimum to move right (higher z) and the minimum value 
of <d(z) decreases. However, because there are small differences between predicted Q(z) in all considered cases then 
verifying the observational test could be difficult. 
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FIG. 5: Confidence levels on the plane (f2 m ,o, f2w,o) minimalized over the rest of parameters for the flat model. The figure 
shows the ellipsoid of the preferred value of n m ,o, fL,o- The results prefer the positive value of Q w ,o, while the negative values 
are allowed (i.e., rotation can exist). 
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FIG. 6: The angular diameter Q for the flat model with rotation for fi m ,o = 0.3 and f2 r ,o = 0.1, 0, —0.02 (top, middle, bottom). 
The minima for these cases are 1.364, 1.605, 1.707, respectively. The rotation causes the minimum to move right (towards to 
higher z) and the minimum value of Q decreases. 
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FIG. 7: The age of the universe to in units of 10 9 j/r for the fiat model with radiation f2 r ,o = 0,0.01 (middle, bottom) and 
rotation f2 w ,o = —0.01 (top). 

V. CONCLUSIONS 

We discuss the problem of universe acceleration driven by global rotation. We demonstrate that the universe 
acceleration increase due to the presence of global rotation effects, although the cosmological constant is still required 
to explain SN la data. In this cosmology, the Friedmann equation is modified by appearance of extra term which 
diminishes with cosmic scale factor as — a~ 4 . Our model suggests limit for rotation fi W) o > —0.033 (at one a level) if 
considered in the Newtonian model, however in the extended model with additional matter which scales like radiation 
(not necessary f2 rj o < 0) we obtain the more safely limit for rotation f2 w .o > —0.01 (at one a level). 

Our limit is weaker than that which can be obtained from BBN (Q Mi q = — 1.23f2 7j o) and CMB (Cl U: o — — 0.41fi 7j o) 
where present value of f2 7j o is estimated as f2 7 .o = 2.48/i~ 2 • 10 -5 |18|. 

We showed that, although the observational constraint from SN la allows only a small contribution from 'dark 
radiation' (however, when in the pure Newtonian model f2 Wj o < a much wider range of negative values of fi^o 
are allowed. We can find the strict analogy between the considered analysis of the observational constraints on the 
global rotation in the model and the search for observational constraints on dark radiation in brane cosmology. The 
corresponding term in brane cosmology scales just like radiation with a constant po or both positive and negative 
Pr,o (Po),o) are possible mathematically. Dark radiation should strongly affect both the Big-Bang nucleosynthesis 
(BBN) and the cosmic microwave background (CMB). Ichiki et al. used such observations to constrain both the 
magnitude and the sign of dark radiation in the case when term p 2 coming from the brane is negligible (it rapidly 
decays as a -8 in the early radiation dominated universe). Therefore, the presence of the term is insignificant during the 
during the later nucleosynthesis. In such an approximation we recover the considered model in which dark radiation 
mimics radiation or rotation. Let us note negative contribution coming from the global rotation presence can reconcile 
the tension between the observed 4 He and D abundance 0|. The application of these results gives also the possible 
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constraints on global rotation term from BBN and from the power spectrum of CMB anisotropies. 

We obtain the limit for Q ufi from BBN as -7.21 • 1CT 5 , while the limit from CMB is -2.41 • 1(T 5 . The present 
extragalactic data suggest uj — 6 • 10 _21 rad/s 0. This gives the strongest limit for fi Wi0 > —5.3 ■ 10 -6 . Therefore, we 
can conclude that the present observational data of SN la give the weaker limit for rotation then obtained by other 
methods. However, let us note that the obtained limitations are constructed in independent manner. 
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